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y  Abstract 

The  material  configuration  considered  is  a  thin  strained  layer  embedded  between  two  relatively 
thick,  unstrained  layers.  The  interfaces  between  the  strained  layer  and  the  unstrained  layers  are 
assumed  to  be  epitaxial,  and  the  strain  in  the  layer  is  assumed  to  arise  from  a  mismatch  in  lattice 
parameter&^tween  the  strained  layerjmate^rial  and  .the  surrounding  material.^  For  purposes  of 
calculation,  the  crystalline  materials  are  assumed  to  have  cubic  symm^^,~^Iented  so  that  the 
[001]  directibfi  is  normal  to  the  interfaces.  Initially,  along  straight  6(^<uslo^tion  exists  on  a ^111^ 
plane,  extending  from  one  otherwise  strain-free  region,  through  the  strained  layer,  and  into  the  other 
strain-free  region.  Thc-purpose-here-is-  t-o  mod^  ^e  process  of  glide  of  this  dislocation^under  the 
action  of  the  stress  associated  with  the  mismatch  strain  in  the  layer.  The  shape  of  the  dislocation 
during  glide  is  approximated  by  a  continuous  piecewise  linear  curve,  and  the  location  of  the  vertices 
of  this  curve  are  viewed  as  generalized  coordinates  of  the  system.  A  corresponding  generalized  force 
is  defined  at  each  node  as  the  variation  in  total  energy  of  the  system  due  to  variation  in  position 
of- that  node.  A  kinetic  law  of  motion,  suitable  for  thermally  activated  glide  of  the  dislocation, 
is  then  cast  into  a  form  involving  the  nodal  velocities  and  the  corresponding  generalized  forces. 
The  governing  equations  have  been  integrated  numerically  for  material  parameters  representative 
of  the  GeSi/Si[001]  strained  layer  systems.  For  layer  thickness  less  than  the  critical  thickness,  the 
dislocation  bows  but  its  total  motion  is  limited.  On  the  other  hand,  when  the  layer  thickness  exceeds 
the  critical  thickness,  glide  reaches  a  stea^  state  after  the  dislocation  in  the  layer  travels  a  distance 
only  several  times  the  layer  thickness.^Preliminary  results  are  also  reported  for  the  situation  when 
the  gliding  dislocation  encounters  an  i  nterface  misfit  dislocation  in  its  path,  the  second  dislocation 
having  arrived  at  the  interface  throug  i  earlier  glide  on  an  intersecting  glide  plane. 
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1.  Introduction 


Due  to  advances  in  the  techniques  of  crystal  growth  over  the  past  decade,  it  is  possible  to  fabricate 
almost  perfect  bi-material  crystals.  If  a  thin  layer  of  material  is  to  be  grown  onto  the  plane  surface 
of  a  second  material  of  the  same  crystal  structure  and  orientation,  atoms  can  be  deposited  in  such 
a  way  as  to  continue  the  crystalline  structure  of  the  substrate  into  the  layer.  Thus,  the  interface 
is  coherent  or  epitaxial,  at  least  in  the  early  stages  of  growth.  Furthermore,  many  layers  can  be 
deposited  in  sequence,  creating  a  material  with  a  periodic  structure  known  as  a  ‘superlattice’.  A 
general  introduction  to  this  technology  and  its  applications  to  electronics  is  given  by  Bean  (1985). 
The  growth  processes  are  sophisticated,  but  are  routinely  performed  by  one  of  several  techniques 
(for  example,  molecular  beam  epitaxy  or  organometallic  vapor  phase  epitaxy)  in  many  laboratories. 


By  proper  selection  of  materials  or  alloys,  it  is  possible  to  match  the  lattice  spacing  of  the 
layer  to  that  of  the  substrate.  Two  examples  are  the  growth  of  AlGaAs  on  GaAs  or  Ino.53Gao.47 As 
on  InP.  In  this  case,  the  layer  grows  without  intrinsic  strain.  These  lattice  matched  structures 
have  been  important  in  the  development  of  electronic  and  opto-electronic  devices  with  desirable 
characteristics,  typically  high  speed  response,  small  size  and  low  power  consumptipn. 


In  the  past  few  years,  the  prospect  of  epitaxial  growth  of  a  layer  on  a  substrate  in  a  situation 
where  the  lattice  parameter  of  the  layer  differs  from  that  of  the  substrate  has  attracted  great 
interest.  A  modest  amount  of  mismatch  is  accommodated  by  an  extensional  strain  in  the  layer,  so 
that  the  crystal  lattice  in  the  layer  remains  perfectly  registered  with  that  of  the  substrate.  This 
process  has  come  to  be  known  as  strained  layer  epitaxy.  Not  all  of  the  effects  of  strain  on  the 
electronic  structure  of  the  material  are  known.  However,  it  appears  that  the  mechanical  strain  can 
be  exploited  to  tailor  the  electronic  energy  band  structure  of  the  material  for  specific  purposes. 
For  example,  strain  can  be  used  to  separate  multiple  modes  at  a  certain  energy  level,  to  convert 
an  indirect  band  gap  material  to  a  direct  band  gap  material,  or  to  adjust  the  size  of  the  energy 
gap  between  the  highest  bound  states  and  the  lowest  conduction  states  of  cliarge  carriers  in  doped 
semiconductor  materials. 


A  difficulty  with  manufacturing  reliable  strained  layer  structures  is  that  the  mismatch  strain 

gives  rise  to  a  driving  force  for  the  nucleation  and  growth  of  defects  in  the  strained  layer  lattice. 

Such  defects  usually  degrade  the  electrical  or  optical  performance  of  the  material,  so  tlie  means 

of  controlling  them  has  become  a  concern  of  central  importance.  A  defect  that  has  attracted 

particular  attention  is  the  threading  dislocation  in  a  strained  layer.  For  a  surface  layer,  tl\is  is  a 
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slip  dislocation  with  the  dislocation  line  extending  across  the  strained  layer,  and  then  continuing 
either  along  the  interface  or  into  the  substrate.  For  an  embedded  layer,  the  line  of  the  threading 
dislocation  may  continue  along  both  interfaces  or  extend  into  adjacent  material  above  and  below 
the  layer.  This  defect  is  important  because  it  can  move  through  the  layer  under  the  action  of 
the  internal  stress,  leaving  behind  segments  of  misfit  dislocation  in  the  interface.  As  the  length  of 
interfacial  dislocation  increases,  the  mismatch  strain  in  the  layer  is  progressively  relaxed. 

The  concept  of  a  critical  layer  thickness  for  the  onset  of  strain  relaxation  has  played  a  central 
role  in  studies  of  epitaxial  layers  (Matthews  et  al,  1970).  There  is  a  simple  explanation  for  this 
phenomenon.  Suppose  that  a  dislocation  extends  from  the  free  surface  of  the  layer  across  its 
thickness,  and  continues  along  the  interface  as  a  misfit  dislocation.  If  the  layer  thickness  is  less 
than  a  certain  value  he  the  work  necessary  to  form  a  length  of  misfit  dislocation  exceeds  the  energy 
released  in  relaxing  the  strain  in  the  layer,  so  the  threading  dislocation  does  not  tend  to  advance 
spontaneously.  On  the  other  hand,  if  the  layer  is  thicker  than  he  the  opposite  is  true  and  there  is  a 
driving  force  tending  to  advance  the  dislocation.  The  discriminating  thickness  is  called  the  critical 
thickness  of  the  layer.  For  a  given  uniform  mismatch  elastic  strain  in  the  layer,  the  threading 
dislocation  will  tend  to  move  so  as  to  reduce  the  length  of  interface  misfit  dislocation  if  the  layer 
thickness  is  below  the  critical  thickness  he  whereas  it  will  tend  to  increase  the  length  of  interface 
misfit  if  the  layer  thickness  exceeds  he 

One  of  the  more  fruitful  extensions  of  this  idea  iias  been  to  estimate  the  spacing  between 
dislocations  in  a  partially  relaxed  layer  at  a  given  initial  mismatch  strain  and  layer  thickness.  This 
extension  is  based  on  the  additional  assumptions  that  the  remaining  elastic  strain  in  the  film  is 
spatially  uniform,  and  that  there  is  a  ready  source  of  dislocations.  However,  dislocation  densities 
observed  experimentally  are  much  lower  than  this  calculation  predicts.  Several  assumptions  in 
the  theory  may  not  be  satisfied  in  practice.  For  example,  the  model  presumes  the  availability  of  a 
supply  of  threading  dislocations  whereas,  in  reality,  dislocations  must  be  nucleated,  a  process  which 
is  typically  more  difficult  to  activate  than  steady  glide  of  a  threading  dislocation  (Fitzgerald  ot  al, 
1989;  Eaglesham  et  al,  1989;  Tuppen  et  al,  1989).  Dislocation  mobility  in  semiconductor  materials 
is  low,  so  that  the  rate  of  elastic  strain  relaxation  can  also  be  limited  by  the  kinetics  of  dislocation 
motion  (Hull  et  al,  1988;  Dodson  and  Tsao,  1987;  Freund  et  al,  1989).  Another  possible  impediment 
to  the  glide  of  threading  dislocations  is  the  interaction  between  dislocations  on  intersecting  glide 
planes.  For  example,  for  the  case  of  350  A  Geo.25Sio.75  layers  on  Si  it  has  been  reported  that  the 
lengths  of  individual  interface  dislocation  lines  did  not  increase  significantly  during  relaxation,  even 
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though  the  number  of  dislocations  increased  (Hull  et  al,  1989).  Furthermore,  the  hesitation  of  a 
threading  dislocation  as  it  encounters  a  misfit  dislocation  in  its  path  is  evident  in  the  real-time  video 
recordings  of  threading  dislocation  motion  in  SiGe/Si  films  during  an  anneal  (Hull  and  Bean,  1990). 
It  was  suggested  that  this  behavior  may  be  due  to  the  interaction  of  dislocations  on  intersecting 
glide  systems. 

The  purpose  here  is  to  examine  the  kinetics  of  glide  of  a  threading  di.slocation  in  an  embedded 
strained  layer.  The  glide  is  assumed  to  be  a  thermally  activated  process  appropriate  for  silicon, 
germanium  and  other  covalent  materials.  Previous  analyses  of  this  process  have  been  based  on 
steady  glide  of  a  dislocation  of  fixed  shape.  In  this  work,  the  complete  shape  of  the  threading 
dislocation  is  calculated  during  both  the  early  transient  and  steady  state  phases  of  its  motion.  The 
three-dimensional  dislocation  line  is  approximated  by  a  series  of  straight  segments,  as  discussed 
in  section  2.  The  shape  of  the  dislocation  is  determined  by  calculating  the  energetic  driving  force 
on  each  segment,  and  integrating  the  resulting  kinetic  equations,  as  discussed  in  sections  3  and  4. 
The  results  for  glide  of  a  single,  isolated  threading  dislocation  in  an  embedded  strained  layer  are 
reported  in  section  5,  and  the  interaction  of  this  dislocation  with  a  pre-existing  dislocation  on  an 
intersecting  glide  plane  is  considered  in  section  6. 

2-.  Description  of  the  Model 

Consider  the  embedded  strained  layer  shown  in  Figure  1.  A  dislocation  is  assumed  to  glide  on  a 
particular  plane  (shaded  in  the  figure).  A  rectangular  xi,a:2i2^3~coordinate  system  is  introduced 
with  orign  in  the  midplane  of  the  layer  as  shown.  The  unit  vector  n  is  the  normal  to  the  glide 
plane  in  the  2:3 -direction.  Suppose  that,  in  the  absence  of  the  gliding  dislocation,  an  equilibrium 
stress  field  cr^  exists  in  the  bimaterial  composite.  This  field  can  be  a  spatially  uniform  isotropic 
membrane  tension  confined  to  the  layer  or  a  nonuniform  stress  distribution  due  to  a  dislocation  on 
an  intersecting  glide  plane.  These  two  situations,  for  the  case  of  a  layer  on  a  single  substrate,  have 
been  recently  considered  in  detail  by  Freund  (1989). 

It  is  assumed  here  that  the  dislocation  tends  to  move  so  as  to  decrease  the  mechanical  energy 
of  the  system.  The  elastic  energy  associated  with  any  particular  dislocation  configuration  is  defined 
as  the  net  work  done  in  introducing  the  fully  formed  dislocation  into  the  stressed  body.  If  is 
the  equilibrium  stress  field  due  to  the  dislocation  in  the  absence  of  any  applied  loading  and  if 
is  the  applied  stress  field,  the  energy  is 

E{S)  =  i  /  b  •  (o-^  .  n)  </5  -1-  /  b  •  (o-^'  •  n)  f/5,  (2.1) 

2  Js  Js 
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where  S  is  the  slipped  area  with  unit  normal  n  and  b  Is  the  Burgers  displacement  vector  of  the 
dislocation.  For  a  glide  dislocation  b  •  n  =  0.  The  boundary  of  the  slipped  area  in  the  glide  plane 
is  denoted  by  B,  and  the  outward  normal  to  B  in  the  plane  is  m.  The  stress  field  cr^  is  sometimes 
called  the  self-stress  of  the  dislocation. 


By  applying  Blin’s  formula  (Ilirth  and  Lothe,  1982)  for  the  energy  of  interaction  between  two 
dislocation  loops,  the  first  integral  in  (2.1)  can  be  expressed  in  terms  of  line  integrals  over  the 
boundary  B  of  the  slipped  area,  resulting  in  the  following  expression  for  the  energy  of  the  system: 

+  £  £(b  X  dll) .  T .  (b  X  dh)  +  (2.2) 

where  R  is  the  separation  distance  between  the  differential  line  elements  dli  and  dl2,  and  T  is  a 
tensor  with  rectangular  components  T,j  =  d^Rldxidxj. 

The  purpose  of  this  section  is  to  describe  a  numerical  method,  based  on  some  implications 
of  (2.2),  for  analyzing  the  change  in  shape  of  a  threading  dislocation.  For  purposes  of  numerical 
simulation,  the  dislocation  configuration  of  Figure  1  is  approximated  by  a  continuous,  piecewise 
linear  configuration.  The  energy  of  the  system  can  then  be  expressed  in  terms  of  the  positions 
of  the  vertices  of  the  piecewise  linear  dislocation.  Suppose  that  this  approximation  consists  of  a 
finite  number  of  straight  dislocation  segments  labeled  Cj,  i  =  0,...,jY  +  1,  each  with  the  same 
Burgers  vector  b.  The  slipped  area  of  this  configuration  is  denoted  by  §  and  its  boundary  by 
B.  The  approjdmate  configuration  is  then  fully  described  by  the  endpoints  of  these  segments. 
Perhaps  the  simplest  description  of  the  configuration  is  illustrated  in  Figure  2.  It  consists  of  jY  -f- 1 
nodes,  N  straight  segments  Ci  of  variable  (but  finite)  length  L,,  and  two  “infinite”  segments  Co 
and  Cn+1-  A  single  generalized  coordinate  Oj  is  cissociated  with  each  node,  and  it  is  a  measure 
of  the  distance  of  the  node  in  the  a:2“direction  from  a  datum  parallel  to  the  xi-axis.  Nodes  are 
constrained  against  motion  in  the  ii-direction.  The  so-called  infinite  segments  are  necessary  ixi 
order  to  properly  model  an  arbitrarily  long  threading  dislocation,  but  they  are  unnecessary  when 
modelling  any  finite  configuration. 


For  the  case  of  a  piecewise  linear  dislocation  configuration,  the  integrals  in  (2,2)  must  be 
evaluated  around  B,  the  boundary  of  the  approximate  configuration.  This  leads  to  the  discretization 
of  the  integrals  in  (2.2) 

AT+l  /V+1 
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N+1  N+1 


E  E  £  £  (*>  X  ■  T  •  (b  X  dk)  +  /  b  ■  (aX  .  „)  dS.  (2.3) 
The  explicit  expressions  for  the  integrals  in  (2.3)  will  be  addressed  in  the  next  section. 


The  variation  of  energy  with  a  local  change  in  the  shape  of  the  dislocation  line  may  be  thought 
of  as  a  force  tending  to  distort  the  dislocation.  The  force  is  defined  as  the  negative  gradient  of  the 
of  the  interaction  energy  with  respect  to  the  position  of  the  dislocation.  Since  the  approximate 
configuration  is  characterized  by  generalized  corner  coordinates  Oj,  it  is  natural  to  introduce  the 
generalized  nodal  force 

^  _  dEjS) 


The  generalized  force  G,-  is  work-conjugate  to  a,  and  acts  in  the  direction  of  the  displacement  o;. 
The  notion  of  a  generalized  force  on  a  dislocation  serves  as  a  convenient  device  for  describing  changes 
in  its  interaction  energy  and  has  been  exploited  in  the  elastic  continuum  theory  of  dislocations  (see 
Dundurs  (1969)).  For  instance,  a  dislocation  configuration  is  said  to  be  in  equilibrium  when  its 
configurational  energy  is  stationary  or,  in  other  words,  when  the  force  on  the  dislocation  vanishes 
pointwise  along  the  dislocation  line.  In  this  investigation,  the  force  is  used  to  calculate  the  rate  of 
change  of  the  corresponding  coordinate  o,-  in  order  to  study  glide  as  a  kinetic  process. 


According  to  Alexander  (1986),  the  process  of  glide  in  covalent  crystals,  such  as  germanium 
and  silicon,  is  a  thermally  activated,  and  the  normal  velocity  of  the  dislocation  line  is  given 
by  a  rate  law  of  the  form 


(2.5) 


where  Vo  is  a  material  constant  with  dimensions  of  speed,  r  is  the  resolved  shear  stress  in  the 
direction  of  the  Burgers  vector  on  the  glide  plane,  Q{t)  is  the  stress  dependent  activation  energy 
for  glide,  k  is  Boltzman’s  constant,  and  T  is  absolute  temperature.  The  notation  =  v  •  m  is 
intended  to  make  clear  the  connection  between  the  dislocation  velocity  v  and  the  normal  m  to  the 
boundary  B  of  S.  A  particular  form  of  the  stress  dependence  of  the  activation  energy  in  covalent 
crystals  has  been  proposed  by  Alexander  (1986)  as 


Q{T)  =  Qo-Eo\og{T/[i),  (2.6) 

where  Qo  and  Eo  are  material  constants.  Values  of  Vp,  Qoi  and  Eo  for  GeSi  alloys  based  on  an 
elementary  mixture  relation  have  been  given  by  Freund  (1989),  and  these  values  will  be  used  in 
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section  5.  At  this  point,  it  should  be  noted  that  the  factor  EojkT  is  of  order  unity  for  the  materials 
being  discussed  here,  and  it  is  assumed  hereafter  that  EolkT  =  1. 

The  form  of  a  rate  law  for  glide  of  the  discretized  dislocation  configuration  in  the  Cj-direction 
at  speed  Vj,  corresponding  to  (2.5),  should  have  the  feature  that  the  nodal  speeds  Vi  depend  on  the 
nodal  forces  G,'  in  the  same  way  that  v  depends  on  r  in  (2.5).  A  local  kinetic  equation  relating  u,- 
to  Gi  and  having  these  properties  is 

Vi  =  (2.7) 


where  the  multiplying  factor  Vo{i)  for  each  node  is  unknown.  The  factors  Vo(i)  are  determined 
by  requiring  that  the  energies  dissipated  through  dislocation  glide  according  to  the  discrete  and 
continuum  descriptions  are  identical,  that  is. 


dl=^ViGu 

i 


(2.8) 


where  B  is  the  boundary  of  5,  and  F{t)  is  the  so-called  Peach-Koehler  force  per  unit  length  on  the 
dislocation  line  in  the  direction  normal  to  the  dislocation  line.  The  Peach-Koehler  force  is  given  by 
F  =  (b  •O’)  X  where  ^  a  unit  tangent  to  and  in  the  direction  of  the  dislocation  (Dundurs  (1969)). 
Thus,  (2.8)  can  be  rewritten  in  the  more  useful  form 


dC=^ViGi. 

i 


(2.0) 


Equation  (2.9)  enforces  the  notion  that  the  generalized  nodal  forces  G,  that  drive  the  glide  of  the 
threading  dislocation  are  the  forces  that  are  work  rate-conjugate  to  the  the  kinematic  rates  i',  with 
respect  to  the  time  rate  of  change  of  the  interaction  energy  E{S). 


To  illustrate  the  ideas  presented  above,  consider  the  dislocation  configuration  shown  in  Figure 
3.  This  simplest  of  piecewise-linear  configurations  consists  of  five  segments,  Co  through  C4,  and 
four  nodes,  gq  through  03,  and  it  is  intended  to  approximate  the  threading  dislocation  of  Figure  1 
gliding  down  the  embedded  layer.  The  segments  Co  and  C4  are  semi-infinite  straight  dislocations 
that  extend  indefinitely  away  from  the  embedded  layer.  The  segments  Ci  and  C3  model  the  misfit 
dislocations  that  are  laid  down  along  the  bi-material  interfaces  as  the  dislocation  segment  C2  is 
driven  down  the  embedded  layer.  The  glide  plane  of  the  dislocation  is  inclined  at  an  angle  a  from 
the  normal  to  the  bi-material  interface.  The  elastic  mismatch  strain  in  the  layer  is  assumed  to  be 
an  isotropic  extensional  strain  Co  in  all  directions  parallel  to  the  interface. 
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To  simplify  the  analysis,  assume  that  the  two  nodes  «o  and  03  remain  fixed  and  that  the  two 
inner  nodes  a\  and  are  constrained  to  move  together  as  one.  The  only  degree  of  freedom  in  this 
problem  is  the  position  o  of  the  center  segment.  The  length  L  of  this  segment  is  the  projected  layer 
thickness  /i/cosa.  It  can  be  shown  that  the  change  in  energy  associated  with  the  self-stress  of  the 


dislocation  due  to  formation  of  the  bow  is 


=  ^[6!  +  [s/FT^-L  -  «  + 


(2.11) 


where  p  is  the  core  cutoff  parameter.  Likewise,  it  can  be  shown  that  the  contribution  to  the  change 
in  energy  due  to  the  stress  Co  induced  by  the  strain  mismatcli  fo  is 


AfJj  =  -2/i7^-^^-^aL6i€oSinacosa. 
(l-i/) 


(2.12) 


The  generalized  nodal  force  G  on  the  dislocation  can  be  found  as  -dAEIda  according  to  (2.4), 


which  leads  to 


G  =  A  [6?  +  -M-1  fl - ]  -  A  +  -S_l  logf — 

27r[^  l-i'jL  L  +  y/W+^i  ^\/)[a 


-j-  \/Z(^  -f  O'] 


(1  -f-  v) 

-b  2u~ - kLbi^o  sin  a  cos  a. 

(l-i/) 


(2.13) 


A  critical  thickness  condition  for  this  approximate  dislocation  configuration  is  obtained  from 
(2.13)  by  requiring  that  G  0  as  a  ->  00.  After  some  elementary  manipulations,  it  is  found  that 
the  critical  layer  thickness  he  for  stability  of  the  threading  dislocation  is  given  by 


Co  =  -rT\ - r; — : - T'l^i  +  (1  “  logf  — 

47r(l -b  1^)710 sin a6i ^  °\2pcosa 

This  result  is  essentially  the  critical  thickness  of  Matthews  ct  al  (1970). 


(2.14) 


In  light  of  (2.14)  the  steady-state  driving  force  Ggt  on  the  dislocation  for  layer  thicknesses 
exceeding  the  critical  thickness  can  be  expressed  as 


+  T^]  [(^)  -  ‘“6 


he 

2p  cos  Q' 


2/)  cos  a /J  ’ 


(2.15) 
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Alternatively,  this  driving  force  can  be  cast  as  an  effective  resolved  shear  stress  (Tsao  and  Dodson, 
1988) 

Gat  f,  fhc\\og{h/2pcosa)\  . 

\  1  ^  ~  T  r  V;  7^ - \\<ToSmacosa.  (2.16) 

oi(/i/cosa)  (  \  h  /  log(/ic/2/Jcosa)  J 


bi{h/i 


(2.16) 


In  this  last  expression,  the  effective  resolved  shear  stress  appears  to  be  linear  ir  CTo.  However,  this 
is  not  the  case  because  the  critical  layer  thickness  he  is  a  function  of  the  layer  stress  through  the 
critical  layer  thickness  criterion  (2.14). 


Using  equations  (2.5)  and  (2.7)  in  conjunction  with  equation  (2.9)  yields 

(2,17) 

.•=0 

where  is  the  normal  velocity  of  the  segment  Ci.  For  this  configuration,  the  only  nonzero 
normal  velocity  is  that  of  segment  C2  and  it  is  given  by  v{a)  =  vi  =  V2i  the  velocity  of  nodes  ci 
and  02 .  Therefore,  (2.17)  simplifies  to 

I  •’(«)«<  =  2V,/.  Jj,  (2.18) 

from  which  it  follows  that 

Vo(2)  =  2  Vo  cos  Of.  (2.19) 


The  velocity  of  the  dislocation  is  therefore 


v{a)  =  2Vc 


G  ] 
jib{h/  cos  a)j 


-Qo/kT 


(2.20) 


where  G  is  given  in  (2.13).  This  velocity  reveals,  among  other  things,  the  rate  at  which  the  threading 
dislocation  lays  down  the  two  misfits  of  length  a.  If  /i  =  he  the  velocity  is  zero.  For  h  >  he,  the 
steady  state  velocity  at  which  the  glide  dislocation  lays  down  the  interfacial  misfits  is  obtained  from 
(2.20)  by  setting  G  =  Gat  in  (2.20).  In  terms  of  the  effective  resolved  shear  stress  derived  earlier, 
the  steady  state  velocity  is  given  by 


(2.21) 


It  is  straightforward  to  extend  this  procedure  to  a  dislocation  line  with  an  arbitrary  number  of 
degrees  of  freedom.  This  issue  is  addressed  below. 
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3.  The  Nodal  Driving  Forces 


In  this  section,  explicit  expressions  for  the  rate  of  change  of  the  elastic  energy  of  the  dislocation 
with  respect  to  variations  in  the  positions  of  the  nodes  are  obtained.  The  basis  for  this  development 
is  found  in  the  expressions  for  the  dislocation  energy  of  two  straight  dislocation  segments  (Hirth 
and  Lothe,  1982).  These  expressions  are  explicit  equations  for  the  first  two  integrals  in  (2.3),  and 
they  are  different  depending  on  whether  the  dislocation  segments  C,-  and  Cm  are  such  that  i  7^  m 
or  i  =  m.  These  results  are  summarized  next. 


With  attention  restricted  to  planar  configurations,  consider  two  dislocation  segments  C,  and 
Cm  which  are  part  of  the  piecewise  straight  dislocation  configuration  of  Figure  2,  oriented  as  shown 
in  Figure  4.  A  local,  oblique  x,  y-coordinate  system  is  generated  from  two  unit  vectors  and 
aligned  with  the  segments  Ci  and  Cm,  respectively.  The  positive  senses  of  these  two  base  vectors 
are  determined  by  requiring  the  unit  vector 


(3.1) 


to  be  directed  in  the  positive  sa-direction.  The  ori^n  coincides  with  the  point  of  intersection  of 
the  lines  x  and  y,  and  0  is  the  angle  between  and  In  this  coordinate  system,  the  segment 
Ci  has  endpoint  coordinates  Xi-i  and  Xi,  and  the  segment  Cm  has  endpoint  coordinates  j/m-i  and 
2/m»  It  is  clear  that  the  x  and  y  axes  lie  on  the  glide  plane  of  the  threading  dislocation,  that  is,  the 
®i>®2“plane;  see  Figure  1. 


It  is  convenient  at  this  point  to  introduce  the  definitions 


Rix,y)  =  \/a;^  +  y*  -  2xy  cos  9, 

(3.2) 

s{x,y)  =  R{x, 7j)  +  y  -  xcos9, 

(3.3) 

t{x,  y)  =  R{x,  y)  +  X  -  y  cos  9. 

(3.4) 

In  terms  of  oblique  coordinates,  the  dislocation  energy  Etm  of  two  straight  dislocation  segments  C, 
and  Cm,  where  i  ^  m,  is  given  by 


where 


I{x,y)  =  xlog 


s{x,y) 

X 


+  ylog 


y 


(3.6) 
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and  I'i^Xatyp)  is  a  special  summation  of  values  of  I{x,y)  at  the  endpoints  of  Ci  and  Cm  defined  as 

•f(®a>  yp)  —  •  ym—1 )  “  1)  J/m)  “  2/m— 1 )  d*  2/m)*  (^•'^) 

The  dislocation  energy  associated  with  a  dislocation  segment,  say  Ci,  interacting  with  itself  (t  =  m) 
is  given  by 

f  [(•-ei)‘  +  l^]i.[log(^)  -l],  (3.8) 

where  p  is  the  core  cutoff  parameter,  and  Li  is  the  segment’s  length. 

The  expressions  (3.5)  and  (3.8)  for  Eim  and  Eu  can  be  used  to  determine  pointwise  the 
dislocation’s  contribution  to  the  energy  E(S)  at  each  node  Oj  by  means  of  the  explicit  summations 
in  (2.3).  Since  each  node  Oj  is  shared  by  two  dislocation  segments  Ci  and  Ci+i,  the  contributions 
to  the  dislocation  energy  at  Uj  can  be  assembled  from  three  sources.  First,  consider  the  interaction 
of  the  segments  Ci  and  C.+i  with  segments  Cm  when  m  <  i,  that  is 

t-i  «-i 

Ejxm  +  .^(t+l)m*  (3*9) 

m=0  m=:0 

Similarly,  the  case  m  >  i  and  the  case  m  =  i,  which  corresponds  to  the  interaction  of  the  two 
segments  sharing  the  node  Oj  must  be  considered.  Finally,  the  contribution  to  the  dislocation  energy 
at  a  node  associated  with  the  self-interaction  of  the  segments  must  be  included.  The  applied-stress 
contribution  to  E{S)  is  obtained  by  numerically  integrating  the  third  integral  in  (2.3)  at  each  node. 

To  determine  the  nodal  driving  force  G,  at  a  given  node  a,,  the  derivative  of  ^(a,)  with  respect 
to  Cj  must  be  formed  according  according  to  (2.4).  This  differentiation  can  be  carried  out  on  each 
individual  contribution  to  E{a,),  such  as  The  driving  force  G,  can  then  be  obtained  from  an 
assembly  operation  analogous  to  that  which  was  described  above  for  the  dislocation  energy.  The 
individual  contributions  to  Gi  are  derived  in  the  Appendix.  Results  are  included  for  nonadjacent 
segments  (A. 11)  and  (A. 14),  adjacent  segments  (A. 17),  segments  interacting  with  semi-infinite 
segments  (A.25),  segments  interacting  with  themselves  (A.26),  and  the  applied  stress  contribution 
(A.29). 


11 


4.  Glide  Kinetics  and  Time  Integration 


In  the  previous  section  expressions  have  been  derived  for  the  nodal  forces  Gi.  These  forces  reflect 
the  tendency  of  the  dislocation  to  change  its  shape  in  order  to  decrease  its  elastic  configurational 
energy.  The  nodal  forces  can  be  used  in  conjunction  with  the  local  kinetic  equation  (2.7)  to  obtain 
the  nodal  velocities  v,-  during  glide  along  the  embedded  layer. 


Consider  once  again  the  piecewise  straight  dislocation  configuration  of  Figure  2.  Restricting 
attention  for  the  time  being  to  the  finite  portion  of  this  configuration,  that  is,  excluding  the  two 
infinite  segments,  the  appropriate  expression  for  each  unknown  coefficient  Vo{i)  at  the  i-th  node 
in  (2.7)  are  derived  first.  To  this  end,  substituting  for  r  from  (2.5)  and  for  Gi  from  (2.7)  in  (2.9) 
results  in 


(4.1) 


where  this  expression  holds  strictly  for  the  case  of  EojkT  —  1.  For  the  approximate  configuration 
B  =  Gi,  (4.1)  can  be  rewritten  as 

JL  rL>  _ 

(4.2) 


where  is  the  normal  velocity  of  segment  of  C,  with  outward  normal  m,,  and  the  arclength  I 
is  taken  as  positive  when  moving  from  node  a,_i  towards  node  a,.  For  a  given  straight  segment 
Cj,  the  normal  velocity  in  terms  of  the  nodal  velocities  and  as  a  function  of  position  along 

the  segment  is 

^  (.1.3) 


where  9i  is  the  angle  between  the  segment  C,-  and  the  horizontal  datum.  Therefore, 


2  ’  [«("*•)(£)] 


(U  -  ^  COS'  +  Vi^iVi  +  vf). 

o 


Then,  solving  equation  (4.2)  for  the  nodal  unknowns  Vo{i)  yields 

ZVoh 


Vo{i)  = 


Li  cos^  6i  +  Li+i  cos^  0i+i 


\<i<N  -  1. 


(4.4) 


(4.5) 


Up  to  this  point,  the  contributions  of  the  two  infinite  segments  Co  and  C'at+i  liave  been 
excluded.  The  results  can  be  extended  to  include  these  contributions  by  using  (4.5)  for  the  cases 
z  =  0  and  i  =  iV,  and  taking  the  limit  Lq  oo  and  Ta'+i  oo  to  obtain  the  trivial  result 

Vo(0)  =  Vo{N)  =  0.  (4.6) 
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Equations  (4.5)  and  (4.6),  in  conjunction  with  the  results  of  section  3  for  the  nodal  forces  G,,  fully 
characterize  the  nodal  velocities  Uj  as  given  by  the  kinetic  equation  (2.7).  It  remains  to  integrate 
these  velocities  in  time  in  order  to  update  the  nodal  coordinates  a; 

At  this  point,  it  is  convenient  to  express  the  driving  force  and  velocities  in  dimensionless  form. 
To  this  end,  introduce  the  non-dimensional  variables  5;,  6,  and  p,  through 

tti  =  aai,  b  =  ab,  =  (4.7) 

where  a  and  7  have  dimensions  of  length  and  stress,  respectively.  Then,  the  elastic  interaction 
energy  of  the  piecewise  straight  dislocation  configuration  which  has  the  following  independent 
variable  dependence 

E  =  E(flo, ...,  a, fl/v!  ^>, /^)  (^•^) 

takes  the  following  convenient  reduction  in  dimensional  dependence: 

E(ao,...,ai,...,a^^;6,/i)  =  Q;^7E(a<„...,n, •,...,  A)  (4-9) 


or, 

E  =  a^jE. 

The  definition  of  E  should  be  clear  in  light  of  (4.9). 


(4.10) 


It  is  now  a  simple  matter  to  express  and  u,  in  terms  of  their  non-dimensional  counterparts 
Gi  and  v,-.  For  the  nodal  forces  it  is  found  that 


Gi  =  a-^Gu 


(4.10) 


where  Gi  =  dEfdci.  As  for  the  nodal  velocities,  letting  1?  =  {Wol pb)ex^{-Qo/kT),  which  has 
dimensions  of  speed,  the  expression  is 

Vi  =  dvi,  (4.11) 


where 


Vi  = 


Gi 


Li  cos*  Oi  4-  co.s2 

In  this  last  expression,  the  lengths  of  the  dislocation  segments  have  been  normalized  by  a. 


(4.12) 


From  among  the  various  ways  in  which  the  nodal  velocities  can  be  integrated  in  tim*',  the 
simplest  numerical  time-integration  scheme  has  been  selected  here,  namely. 


=  a’y  +  Atvf  l<i<N-l, 


(4.13) 
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where  the  superscripts  n  and  n  +  1  identify  values  at  consecutive  time  steps.  The  non-dimensional 
time  increment  At  in  (4.13)  is  equal  to  (i?/a)At.  As  with  all  forward  or  explicit  Euler  time 
integration  schemes,  some  sort  of  automatic  time-step  control  is  desirable  in  order  to  obtain  accurate 
results  from  (4.13).  A  suitable  local  error  associated  with  the  integration  at  each  node  a;  is 


,n  iA..vuri 

—  i-nl  1 


(4.14) 


where  Vv”  =  v"  -  that  is,  the  first  backward  difierence  of  u"  (Press  et  al,  1986).  Then  the 
automatic  time-step  control  can  be  achieved  by  controlling  a  global  error  associated  with  the  entire 
configuration,  say 

e”  =  max  lef  1.  (4.15) 


5.  Simulation  of  a  Single  Threading  Dislocation 

In  this  section,  the  motion  of  a  single  threading  dislocation  in  an  embedded  strained  epitaxial  layer 
is  simulated.  The  elastic  mismatch  strain  in  the  layer  is  an  isotropic  extensional  strain  Co  the 
plane  parallel  to  the  interface.  For  definiteness,  the  layer  is  assumed  to  be  the  alloy  Ge^Sii-*,  where 
X  is  the  fraction  of  lattice  sites  occupied  by  germanium  atoms,  and  1  -  x  is  the  remaining  fraction 
of  sites  occupied  by  silicon  atoms.  The  numerical  simulations  are  restricted  to  the  Ge.25Si.75  alloy. 
The  choice  of  material  is  motivated  by  experimental  observations  reported  by  Hull  et  al.  (19S8)  for 
values  of  x  between  0.25  and  0.3.  The  surrounding  material  is  Si  and  the  materials  are  oriented  so 
that  the  normal  to  the  interfaces  is  the  [001]  crystallographic  direction.  The  threading  dislocation 
glides  on  a  {111}  plane  which  is  inclined  at  an  angle  a  =  sin“\l/>/3)  with  respect  to  the  interface 
normal.  In  terms  of  xi,X2~coordinates,  the  components  of  the  Burgers  vector  of  the  dislocation 
are  61  =  a/36/2  and  62  =  -6/2.  The  cutoff  parameter  p  has  been  set  at  6/8.  As  the  threading 
dislocation  glides  down  the  layer,  it  leaves  behind  a  pair  of  straight  60°  interface  misfit  dislocations, 
one  on  each  interface. 

Consider  the  behavior  of  a  dislocation  threading  the  layer.  It  is  assumed  that  the  layer  thickness 
exceeds  the  critical  thickness,  so  that  the  dislocation  is  unstable  in  the  sense  of  elastic  equilibrium 
theory,  and  the  threading  segment  glides  indc.'.aitcly  far  along  the  layer.  Two  basic  questions 
that  may  be  asked  are:  What  is  the  magnitude  of  the  driving  force  on  the  dislocation  during  the 
transient  glide  process,  and  what  is  the  rate  of  misfit  formation?  Tliis  section  addresses  these  two 
issues.  In  addition,  the  shape  of  tlie  threading  dislocation  at  various  stages  of  the  growth  process 
is  calculated. 
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Following  a  suggestion  by  Freund  (1989),  kinetic  data  for  Ge~Sii_i  films  are  estimated  from 
available  data  on  germanium  and  silicon  by  applying  the  rule  of  mixtures.  Therefore  parameters 
such  as  Eo,  Qo,  and  Vo  which  were  introduced  in  section  2  are  determined  from 

Eo-xEg-^{l-x)Ei,  Qo  =  xQg  +  (l-a:)Qj,  Vo  =  xVg-V{l- x)Vs.  (5.1) 

Based  on  data  tabulated  by  Alexander  (1986),  Eg  =  0.064  eV,  Q,  =  1.6961^,  and  V,  =  1.36  x 
10®  m/s  for  silicon  and  Eg  =  0.195  eV,  Qg  =  0.307 cV,  and  =  1.145  X  10®  m/s  for  germanium. 
Therefore,  the  kinetic  data  used  for  Ge.25Si.75  is  Eo  =  9.68  X  10“^  eV,  Qo  =  1.34  eV,  and  Vo  = 
3.88  X  10®  m/s.  For  this  particular  value  of  Eoi  EojkT  varies  from  1.29  to  0.88  over  a  temperature 
range  of  600-1000C,  therefore  the  assumption  that  EojkT  ~  1  is  a  reasonable  one  over  temperature 
ranges  where  experimental  observations  have  been  made. 

Two  case  geometries  have  been  considered,  namely,  h  =  1.1/ic  and  h  =  2/ic.  In  both  cases  the 
mismatch  strain  £<>  =  0.02.  The  simulations  are  carried  out  on  the  glide  plane  of  the  threading 
dislocation,  the  xi,X2-plane,  where  the  uniform  isotropic  stress  ao  =  2/i€o(l  +  i')/(l  -  acting 
in  the  plane  of  the  layer,  induces  a  resolved  shear  stress  acting  in  the  direction  of  the  Burgers 
vector  of  the  dislocation.  This  resolved  traction  is  confined  to  the  projected  layer  thickness  on  the 
glide  plane,  the  thickness  of  which  is  hp  =  h/cosa.  For  h  =  l.lhg*  tbe  normalized  projected  layer 
thickness  hpjh  is  37.49.  For  h  ~  2/ic,  the  normalized  projected  layer  thickness  is  68.17.  Outside 
the  layer,  a  Peierls  stress  of  magnitude  1  x  10~®/t  is  applied  to  simulate  the  lattice  resistance  to 
the  motion  of  the  dislocation.  The  ratio  of  the  Peierls  stress  to  the  resolved  stress  in  the  layer  is 
1:35  for  u  =  0.3. 

The  piecewise  linear  approximation  of  the  threading  dislocation  is  made  up  of  three  sections. 
The  middle  section  -  the  layer  section  -  consists  of  an  odd  number  of  evenly  spaced  nodes,  the 
first  node  being  on  one  of  the  interfaces  and  the  last  node  being  on  the  other  interface.  The  nodal 
spacing  is  measured  along  the  xi-axis  and  is  approximately  equal  to  96  in  both  cases.  The  other  two 
sections  -  the  “wing”  sections  -  model  the  two  sections  of  the  dislocation  line  that  lie  outside  the' 
projected  layer  thickness.  The  ratio  of  the  projected  layer  tliickncss  to  the  length  of  the  dislocation 
wings  is  approximately  1:10  in  both  cases.  The  outer  node  of  a  wing  section  lies  on  the  xi-axis 
and  is  attached  to  one  of  the  semi- in  finite  elements  discussed  in  the  Appendix.  The  semi-infinite 
elements  run  along  the  xj-axis  and  are  not  allowed  to  move  during  the  kinetic  simulation  -  in  effect, 
the  threading  dislocation  is  pinned  at  two  nodes  which  are  relatively  far  from  the  interfaces. 
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At  the  start  of  the  simulation  the  dislocation  is  straight  and  parallel  to  the  xi-axis.  As  the 
dislocation  deforms  and  glides  down  the  layer  due  to  the  elastic  mismatch,  the  dislocation  segments 
grow  in  size  significantly,  particularly  the  element  on  each  wing  closest  to  the  layer.  For  this  reason, 
an  automatic  remeshing  scheme  has  been  devised  which  subdivides  these  two  dislocation  segments 
based  on  their  current  lengths.  The  remeshing  stops  once  the  distance  of  the  closest  wing  node 
to  the  interface  reaches  a  critical  minimum  distance  along  the  ari-axis  and  the  dislocation  segment 
is  nearly  parallel  to  the  interface.  Due  to  the  single  degrce-of-frcedom  formulation  of  the  model, 
a  dislocation  segment  that  is  parallel  to  the  a;2-axis  would  be  inadequate  and  would  be  unable  to 
properly  characterize  the  interaction  of  its  two  end-nodes. 

Figure  5  shows  the  motion  of  the  threading  dislocation  on  a  layer  with  thickness  h  =  1.1/ic  at 
an  early  non-steady  stage  of  its  growth  history.  The  figure  shows  several  configurations  starting  at 
time  t  =  5  X  10^  until  time  <  =  1.7  X  10®  at  regular  time  intervals  At  =  1.5  X  10'*.  As  was  shown  in 
section  4,  the  time  scale  used  in  the  simulations  is  non-dimensional  and  real  time  can  be  obtained 
by  multiplying  the  normalized  time  by  (a/3Vo)exp{Qo/kT).  For  a  temperature  T  =  700A'  the 
dimensional  factor  is  1.52  X  10"®s,  for  T  =  900K  it  is  1.09  X  10~®s,  and  for  T  =  1200K  it  is 
1.45  X  10~*°s.  The  higher  the  temperature,  the  greater  is  the  velocity  of  the  dislocation.  The 
asymmetry  in  the  dislocation’s  shape,  which  has  been  observed  by  Hull  and  Bean  (1990),  for 
example,  is  due  to  the  mixed  nature  of  the  Burgers  vector.  If  the  components  of  the  Burgers  vector 
were  bi/b  =  1  and  62/6  =  0,  then  the  dislocation  would  be  symmetrical  about  the  X2-a.xis. 

Figure  6  shows  the  long-term,  steady-state  motion  of  the  dislocation.  Here,  the  time  interval  At 
between  configurations  is  6x  10“*.  The  first  configuration  shown  in  the  figure  is  at  <  =  5.5x  10'*.  The 
length  of  the  misfit  at  t  =  4.75  X  10®,  the  last  frame  in  the  figure,  is  approximately  fifteen  times  the 
projected  layer  thickness.  It  is  apparent  from  the  figure  that  the  velocity  of  the  dislocation  is  nearly 
constant.  Figure  7  shows  non-dimensional  values  of  the  nodal  driving  force  G,  for  various  nodes  in 
the  layer  as  a  function  of  time.  It  is  perhaps  surprising  that  the  forces  acting  on  the  nodes  within 
the  layer  are  different  from  each  other,  even  though  the  nodes  are  all  moving  at  essentially  the  same 
speed.  Although  all  nodcil  forces  were  normalized  in  the  same  way  in  the  kinetic  equation  (2.7), 
the  coefficient  Vo{i)  at  any  particular  node  depends  on  the  configuration  of  the  segments  meeting 
at  that  node  (cf.  (4.5)).  It  is  possible  that  some  other  scheme  for  normalizing  G,  in  (2.7)  would 
result  in  the  forces  on  nodes  being  equal  when  the  nodal  velocities  are  equal,  but  such  a  scheme  has 
not  yet  been  identified.  Figure  8  shows  the  same  information  as  a  function  of  nodal  position.  The 
slight  discontinuity  in  the  figures  (near  t  =  5  in  Figure  7,  say)  is  due  to  the  remeshing  procediue 
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described  above.  To  obtain  physical  values  from  the  non-dimensional  force  one  must  multiply  the 
simulation  data  by  0^7,  as  shown  in  section  4.  Assuming  that  7  =  7.08  X  10^dynes///m^  (silicon) 
and  a  Burgers  displacement  of  4A,  this  factor  is  1.13  x  10"^dynes.  As  these  two  figures  show, 
initially  the  driving  force  is  quite  high,  particularly  on  the  center  and  inner  nodes.  With  time,  the 
force  sharply  drops  and  reaches  a  slowly  decreasing  plateau.  In  light  of  the  discussion  at  the  end 
of  section  2,  the  driving  force  will  reach  an  asymptotic  steady-state  value  as  i  00.  For  practical 
purposes,  though,  the  simulations  reveal  that  this  force  becomes  nearly  constant  after  the  misfit 
length  is  a  few  times  the  projected  thickness  length.  The  rate  of  misfit  formation,  which  is  given 
by  the  velocity  of  a  node  at  the  interface,  is  about  3.614  exp(-Qo/^'T)  when  the  driving  force  has 
reached  a  steady  level. 

Figure  9  shows  the  early  stages  of  growth  of  a  threading  dislocation  gliding  down  a  layer  with  a 
thickness  h  that  is  2  times  the  critical  layer  thickness  he.  Five  configurations  are  shown,  starting  at 
f  =  i  X  10'*  and  spaced  at  regular  time  intervals  At  =  1  x  10'*.  The  long  term,  steady-state  motion 
of  this  dislocation  is  shown  in  Figure  10,  where  the  first  configuration  is  at  t  =  7.5  X  10'*  and  the 
time  interval  At  between  configurations  is  2.5  x  10'*.  The  length  of  the  misfit  at  t  =  2  x  10®,  the  last 
frame  in  the  figure,  is  approximately  twenty  times  the  projected  layer  thickness.  The  final  time  in 
Figure  10  is  less  than  half  of  that  in  Figure  6.  It  is  apparent  that  the  dislocation  shown  in  Figure 
10  is  gliding  at  a  greater  speed  than  that  of  the  dislocation  shown  in  Figure  6.  This  observation 
is  consistent  with  equation  (2.21)  which  reveals  that  the  steady-state  velocity  Vjt  increases  with 
increasing  driving  force  G.  Figure  11  shows  non-dimensional  values  of  the  nodal  driving  force  G, 
for  some  of  the  layer  nodes  as  a  function  of  time,  and  Figure  12  shows  the  same  information  as 
function  of  position.  The  driving  force  on  the  center  node  is  about  five  times  larger  than  it  is  in 
Figure  8.  Furthermore,  the  steady-state  rate  of  misfit  formation  is  about  22.Wo^x.p{-QolkT), 
which  is  close  to  seven  times  faster  than  before. 
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6.  Simulation  of  Interactions  on  Intersecting  Glide  Planes 


Up  to  now,  the  numerical  simulations  have  been  restricted  to  the  motion  of  a  single  threading 
dislocation  through  an  embedded  layer  which  is  uniformly  strained.  Of  interest  are  cases  were  a 
non-uniform  strain  distribution  exists  on  the  layer  due  to  additional  traction  on  the  glide  plane. 
The  non-uniform  strain  distribution  changes  the  driving  force  acting  on  the  threading  dislocation. 
In  this  section,  one  such  case  is  considered,  namely,  the  blocking  of  a  threading  dislocation  by 
an  interface  misfit  dislocation  lying  on  an  intersecting  glide  plane.  The  specific  configuration  to 
be  considered  is  shown  in  Figure  13.  For  the  remainder  of  this  section,  the  coordinate  axis  labels 
adopted  previously  will  be  dropped  in  favor  of  those  shown  in  the  figure.  The  two  glide  planes  shown 
are  intersecting  {lll}-planes.  These  planes  are  identified  on  the  pyramids  of  close-packed  {111}- 
planes  in  a  cubic  material  in  Figure  14.  As  in  section  5,  the  embedded  layer  and  the  surrounding 
material  are  assumed  to  be  Geo.25Sio.75  and  Si,  respectively,  and  the  normal  to  the  interfaces  is 
the  [001]  crystallographic  direction.  With  respect  to  the  base  vectors  of  the  a‘i,X2,a;3-coordinate 
system,  the  normal  to  the  xj > ®3-plane is  the  vector  {>/273,0,  \/TjZ)  and  the  normal  to  the  xi,X3- 
plane  is  the  vector  {0,  \/2/3,  \/l/3}. 

As  shown  in  Figure  13,  the  threading  dislocation  glides  on  the  X2,X3-plane  in  the  positive 
X2-direction,  and  it  is  assumed  that  its  Burgers  vector  lies  in  this  plane.  A  misfit  dislocation, 
modelled  by  an  infinite  straight  dislocation,  lies  along  the  xj-axis;  its  Burgers  vector  is  assumed  to 
lie  in  the  xi,X3 -plane.  Image  effects  associated  with  the  interfaces  will  be  ignored.  As  discussed 
by  Freund  (1989),  there  are  several  possible  combinations  of  Burgers  vectors,  but  only  cases  were 
both  dislocations  relieve  elastic  strain  in  the  layer  are  of  interest.  The  possible  directions  of  Burgers 
vectors  for  the  two  dislocations  are  represented  by  the  edges  of  the  two  shaded  planes  in  Figure  11. 
These  correspond  to  eight  combinations  of  Burgers  vectors  with  components  in  the  X3-diroction. 
These  combinations  are  reproduced  below  in  Table  1  as  they  were  listed  by  Freund  (1989). 

Tabic  1 
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The  magnitudes  of  the  two  Burgers  vectors  are  assumed  to  be  the  same.  Moreover,  the  components 
of  the  vectors  are  referred  to  unit  base  vectors  of  the  a:i,a;2>®3-coordinate  system. 


The  stress  field  arising  from  the  misfit  dislocation  along  the  xi-axis  sets  up  an  additional 
contribution  to  the  applied  stress  field  in  (2.1).  Denote  this  additional  stress  field  by  for  the 
present  discussion.  This  field  is  given  by  the  stress  field  of  a  straight  dislocation  in  an  unbounded 
body.  Hence,  the  stress  field  corresponding  to  the  edge  component  of  the  misfit  dislocation  depends 
only  on  position  in  the  a:2j®3-plane,  and  it  can  be  computed  from  the  Airy  stress  function 

log  ^xi  +  (»3  +  /./2)= 


X(a:2,a;3) 


47r(l  —  I/) 


(6.1) 


according  to 
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The  stress  field  associated  with  the  screw  component  of  the  Burgers  vector  of  the  straight  dislocation 
is 
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The  above  stress  distribution  gives  rise  to  a  resolved  shear  stress  r  =  [b  •  (<r^  •  n)]/h  in  the 
direction  of  the  Burgers  vector  on  the  X2,a:3-plane.  This  stress  can  increase  or  decrease  the  mag¬ 
nitude  of  the  nodal  driving  force  <j,  on  the  dislocation.  The  interaction  nodal  driving  force  Gf 
is  estimated  by  the  method  of  the  Appendix  A.5.  Instead  of  the  linear  quadrature  rule  used  in 
that  section  to  integrate  the  contribution  of  the  misfit  strain,  a  suitable  quadrature  rule,  that  takes 
into  account  the  fact  that  the  stress  field  cr^  decreases  as  one  over  the  radial  distance  to  the  misfit 
dislocation,  is  used. 
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The  two  layer  geometries  studied  in  section  5,  namely,  h  =  l.lhc  and  h  =  2hc,  have  been 
reconsidered,  and  once  again  the  mismatch  strain  £„  in  the  embedded  layer  is  0.02.  The  numerical 
simulations  are  carried  out  for  the  Burgers  vector  combination  (i)  as  defined  in  Table  1.  For  this 
combination,  the  non-uniform  strain  distribution  encountered  by  the  advancing  glide  dislocation 
is  such  that  its  driving  force  is  reduced  as  it  approaches  the  misfit  dislocation  from  the  negative 
a:2-direction;  see  Freund  (1989, 1990)  for  more  details. 

At  the  beginning  of  the  simulation,  the  threading  dislocation  starts  parallel  to  the  Xg-axis  at 
X2  =  -/,  where  I  is  some  suitable  distance  from  the  orign.  The  straight  misfit  is  always  positioned 
at  X2  =  0  and  X3  =  hflcosa.  For  the  first  case,  h  =  l.lhc,  simulations  where  done  for  values  of 
/  =  150b,  3006,  and  6006.  Figure  15  shows  a  sequence  of  configurations  leading  to  the  blocking  of 
the  threading  dislocation  by  a  misfit.  The  first  configuration  shown  is  at  time  f  =  2.5  x  10'*  and  the 
last  one,  at  t  =  6.5  X  lO"*.  The  time  interval  between  configurations  is  At  =  5  x  10^.  The  threading 
dislocation  was  started  at  X2  =  -1506.  Figure  16  shows  non-dimensional  values  of  the  nodal 
driving  force  (?<  for  the  center  layer  node  and  the  interfacial  node  as  a  function  of  nodal  position 
along  the  X2-axis.  The  dotted  curves  corresponding  to  values  of  Gi  in  the  absence  of  a  misfit  are 
shown  for  comparison.  The  figure  shows  that  the  driving  force  on  the  threading  dislocation  remains 
unaflfected  until  the  dislocation  is  within  a  radius  of  256  from  the  misfit.  At  this  point,  the  near  field 
of  the  straight  dislocation  becomes  dominant  and  decreases  the  driving  force  quite  severely.  This 
slows  down  the  gliding  dislocation  which  ultimately  arrests.  The  numerical  simulation  presented 
here  is  unsuitable  for  predicting  configurational  issues  after  the  intcrfacial  segment  of  the  threading 
dislocation  intersects  the  misfit.  At  this  point,  a  dislocation  reaction  would  take  place,  resulting  in 
the  formation  of  a  jog  or  a  kink.  A  realistic  simulation  would  have  to  take  into  account  the  newly 
created  dislocation  segment. 

In  order  to  illustrate  the  slow-down  of  the  gliding  dislocation  due  to  the  presence  of  a  misfit. 
Figure  17  shows  a  superposition  of  sequences  of  configurations.  The  solid  curves  correspond  to 
a  gliding  dislocation  being  blocked,  the  dotted  curves  correspond  to  a  dislocation  gliding  in  the 
absence  of  a  misfit.  The  initial  configurations  shown  are  at  time  /'  =  1  X  10®  and  the  last  ones,  at 
t  =  1.5  X  10®.  The  time  interval  between  configurations  is  Ai  =  1  x  10'*.  The  threading  dislocation 
was  started  at  X2  =  -3006.  It  is  apparent  from  the  figure  that  the  slow-down  is  significant.  Figure 
18  shows  non-dimensional  values  of  the  nodal  driving  force  G,  for  the  center  layer  node  and  the 
interfacial  node  as  a  function  of  nodal  position  along  the  X2-axls.  Finally,  Figure  19  shows  non- 
dimensional  values  of  the  nodal  driving  force  G,  for  the  center  layer  node  and  the  interfacial  node 
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as  a  function  of  nodal  position  along  the  X2-axis,  for  the  case  where  the  threading  dislocation  is 
started  at  X2  =  -6006.  Here,  the  magnitude  of  the  driving  force  on  the  gliding  dislocation  is  rather 
low  as  the  dislocation  approaches  the  misfit  at  X2  =  0,  so  only  a  small  change  in  G  is  required  for 
the  dislocation  to  arrest. 

Results  for  the  second  layer  geometry,  h  =  2/ic,  are  very  similar  to  those  illustrated  above  for 
h  =  l.lhc.  Figure  20  shows  a  sequence  of  configurations  leading  to  the  blocking  of  the  threading 
dislocation  by  a  misfit.  The  threading  dislocation  was  started  at  X2  =  -3006.  As  discussed 
in  section  5,  for  this  configuration  the  driving  force  is  much  higher  and  consequently  it  is  more 
difficult  to  arrest  nodes  in  the  middle  of  the  layer.  Figure  21  shows  non-dimensional  values  of  the 
nodal  driving  force  Gi  for  the  mid-layer  node,  an  intermediate  node,  and  the  interfacial  node  as  a 
function  of  nodal  position  along  the  X2-axis.  In  contrast  with  the  case  of  /i  =  l.lhc,  tbe  middle 
node  is  barely  affected  by  the  stress  field  of  the  straight  misfit.  Consequently,  the  slow-down  is  less 
dramatic. 
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Appendix:  Expressions  for  Nodal  Driving  Forces 


A.l  Nonadjacent  Interacting  Segments 

Consider  a  segment  with  end  nodes  flj  and  Oj+i  interacting  with  a  second  disjoint  segment, 
say  Cm+ii  with  end  nodes  and  Om+i-  Figure  22  shows  the  geometry  for  this  case.  The  nodal 
coordinates  o^,  Cj+i,  a^j  a^nd  am+i  are  measured  from  a  common  horizontal  datum  as  was  described 
in  section  2.  The  horizontal  distance  between  the  ends  of  the  segment  Cj+i  is  d,+i;  for  the  segment 
Cm+iy  this  distance  is  dm+i*  From  the  geometry  it  is  found  that  the  nodal  coordinates  in  the 
oblique  a:,j/,^-system  are 

®t+l  ~  "b  (■''^•2) 

y^  =  {D-  X)  (.1.3) 


2/m+l  —  Vm  "b 


where 


~  ^m)dt+l  (flm  ~  Q»')di4.;c/,n+i 
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Li+\  =  yj (a, -4.1  -  a,)2  + 

Lm+l  ~  ^ (®m+l  —  ^m)"  "b  d^n+l> 
and  D  is  the  horizontal  distance  between  nodes  a,-  and  a,„. 


In  light  of  (3.5),  the  interaction  energy  between  the  two  segments  is 


M  1 


.7i+l'/m+l 


•^(®a)  I//?)} 


where 


fr  =  drbi+{ar  -  ar-l)b2, 

9r  ~  (*^*r  ~  ®r— l)^l  ~  (Irbl-  ( 

It  follows  that  the  rate  of  change  of  the  dislocation  energy  with  respect  to  nodal  position  is 

^  IP  ^  f  I.  f  t  biflm+l  rf  1,  ^  f  e  I  9i+\9m  +  l 


(.1.9) 

(.1.10) 


(.1.11) 
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In  the  above  expression,  dl{xa,yp)ldai  can  be  written  as  +  *^0),  where  the  functions 

Ha  and  Ja  are  defined  as 


dXa  ,  Saa  ,  ^2/ot  ^aa 

—  log - +  —  log  - - , 

Udx  Udx  t{3— Of)or 

Xg  {  dSga  dSa(3-a)  } 

SaaSg(3-a)  I  da;  dtti  } 

,  Va  f  dt  aa  ,  j  9^(3  — a)a  1 

““  dai  / 


(/1.12) 


(A.13) 


In  these  last  two  expressions,  Sa(3_a)  denotes  the  function  s(a;,y),  as  defined  in  (3.3),  evaluated 
aA  X  =  Xg  and  y  =  3/3-0.  It  has  no  connection  with  the  notation  introduced  in  (3.7)  for  I{x,y). 
Derivatives  such  as  dsgajdai  are  readily  obtained  from  expressions  displayed  in  this  section  and 
the  details  are  omitted. 


Another  instance  of  interaction  of  nonadjacent  segments  consists  of  a  segment  with  end  nodes 
o<_i  and  Oi  interacting  with  a  second  segment  with  end  nodes  and  am+i .  Figure  23  shows  the 
geometry  for  this  case.  The  previous  analysis  holds  for  this  case  as  well  if  a,  is  replaced  by  Oj-i, 
fli+i  by  a,-,  and  so  on,  and  corresponding  algebraic  signs  are  changed.  It  follows  that  the  rate  of 
change  of  interaction  energy  for  this  case  is 
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A.2  Adjacent  Interacting  Segments 


Next,  the  interaction  of  two  adjacent  dislocation  segments  which  share  the  node  a,  is  considered. 
Figure  24  illustrates  this  situation  where  the  nodal  coordinates  have  been  chosen  as  x,  =  0,  x,+i  = 
Li+i,  Vi-i  -  -Li,  and  yi  =  0.  The  lengths  X.+i  and  Li  of  the  x,x,+i  and  yi-iVi  segments  are 
given  by 

T,+i  =  (A. 15) 

Li  =  {ai  -  tti-iY  +  d} ,  (A.16) 


An  analysis  analogous  to  that  presented  in  section  A.l  reveals  that  the  rate  of  change  of  dislocation 
energy  for  adjacent  segments  is 
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(A.17) 


A.3  Segment  Interacting  with  Infinite  Segments 

The  case  of  a  segment  interacting  with  two  semi-infinite  segments  requires  special  attention.  As 
was  mentioned  in  section  2,  the  semi-infinite  segments  are  necessary  to  model  arbitrarily  long 
dislocations,  but  are  not  needed  to  model  finite  dislocation  loops.  Figure  25  shows  a  segment  C,+i 
with  end  nodes  a<  and  Cj+i  and  two  infinite  segments  y'ly't  and  y'/.yV •  figure,  the  length  of  an 

infinite  segment  is  shown  to  be  m,  but  the  limiting  case  m  -*  oo  is  of  primary  interest.  .\  difficulty 
associated  with  determining  the  energy  of  interaction  of  the  dislocation  configuration  shown  in 
Figure  25  is  the  fact  that  this  energy  is  unbounded  as  m  — ‘  cx).  .\  solution  to  this  problem  is  to 
consider  the  energy  increase  in  going  from  a  reference  configuration  consisting  of  segments  y[y'2 
and  y'fy'ft  the  dummy  segment  y,y,+i,  to  one  consisting  of  segments  y'^yt,  x,x,+i,  and  y'/i/o- 
This  is  the  same  procedure  that  was  followed  in  the  exjimple  at  the  end  of  section  2  to  obtain  tlie 
energy  associated  with  the  formation  of  a  kink.  The  energy  change  is 

AjEin/  =  E{y'iy2 :  liajj+i)  +  E(xiXi+i  :  y'ly")  -  E{y\y2  :  ViVi+i)  -  Eiy'iy"  ;  yiVi+i).  (-4.18) 

This  energy  change  is  finite  in  the  limit  as  m  — »  oo  because  the  singularities  associated  with  the 
two  semi-infinite  dislocation  segments  cancel  each  other  out.  Moreover,  the  final  expression  for  G, 
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is  unaffected  by  the  choice  of  AJSi„/  over  that  of  E,nj  because  the  reference  configuration  does  not 
depend  on  the  nodal  coordinate  a,-. 


The  first  term  E{y[y2  :  x.x.+i)  +  :  y'^y”)  is  given  by 

E{y[y2  :  a:,-a:,+i)  +  jE(a:,-x,+i  :  y'^y'^)  =  —  cose  l{xa,yp), 


(4.19) 


where  /(xq,  j//?)  can  be  decomposed  additively  into  a  non-singular  term  /o,tm  and  a  singular  term 
The  implication  is  that  the  former  term  remains  bounded  in  the  limit  as  m  oo,  whereas 
the  latter  term  does  not.  These  two  contributions  can  be  written  explicitly  as 

Io{xa,yi3)  =  a:,-  log  ^  +  Xi+i  log  +  y"  log  +  y'^  log  (4.20) 

and 

Imixa, yp)  =  Xi log ^  -1-  Xi+i  log  +  y[  log  +  ?/(' (4.21) 

^i2  ^(t+l)l  ^'2 

In  equations  (A.20)  and  (A.21)  some  abbreviated  notation  has  been  adopted.  For  instance,  s',^  = 
s{xi,y[)  where  s{x,y)  has  been  defined  in  (A.3). 

The  remaining  terms  in  (A.  18)  are 


EWiy2  :  yiVi+i)  +  E{y'ly'2  :  ym+i)  =  ■j-I{xa,yp)y 


(4.22) 


where  it  can  be  shown  that 

I{xa,yp)  rf.-nlog  +  di+i){D"  -  di+i)  ^  ^D"-di+i 

D'\og-^J'-^>-  +  2di+i. 

The  definitions  of  D" ,  D',  and  d,+i  should  be  clear  from  Figure  25. 

Combining  all  the  terms  in  (A. 18)  and  taking  the  limit  m  -*  oo  yields 

*  T  /  \  ,  2(/,+i  ,  sin 01  . 

AEinf  =  -^COS^I /o(Xa,i//j) -t- -  2.C,+  1  logX,+i  -  (.pgg 

^|(d.+,  -  I»")log(D"-d,-+i)-H?'log;5'-l- 

idi+,-D')\ogD'  +  di+xy 


(4.23) 


(4.24) 
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The  above  result  is  a  generalization  of  equation  (2.11).  Finally,  the  rate  of  change  of  the  dislocation 
energy  with  respect  to  nodal  position  is  found  to  be 
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(A.25) 


In  addition  to  the  interaction  of  a  dislocation  segment  with  two  infinite  segments,  there  are  two 
other  special  interactions  of  nonadjacent  segments  which  must  be  considered  in  order  to  successfully 
implement  the  numerical  scheme.  Both  Ccises  come  about  when  the  coordinate  system  described  in 
the  introduction  of  this  section  breaks  down.  The  first  involves  colinear  segments,  either  adjacent 
or  nonadjacent.  The  second  case  involves  the  interaction  of  parallel  nonadjacent  segments.  Suitable 
expressions  can  be  derived  as  limiting  cases  of  the  equations  derived  in  sections  A.l  and  A.2.  The 
case  of  parallel  segments  is  discussed  by  Birth  and  Lothe  (1982). 


A.4  Self-Interacting  Segment 


Finally,  to  conclude  the  derivation  of  the  contribution  to  the  nodal  forces  associated  with  the  dislo¬ 
cation  self-stress  tr^,  the  dislocation  energy  associated  with  self-interacting  segments  is  considered. 
Once  again,  attention  is  focussed  on  the  geometry  of  Figure  21  which  shows  that  the  node  a,  con¬ 
tributes  to  the  dislocation  energy  E„  of  two  dislocation  segments.  Making  use  of  equation  (3.8) 
and  some  of  the  results  of  section  A.2,  it  is  found  that  the  rate  of  change  of  energy  with  respect  to 
nodal  position  is 
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A. 5  Applied  Stress  Contribution 


The  elastic  interaction  energy  aissociated  with  the  work  done  during  glide  by  the  applied  stresses 
acting  on  the  glide  plane  was  first  introduced  in  equation  (2.1).  This  energy  can  be  written  as 

Ea{S)  =  j  brdS,  {A.27) 

where  the  subscript  A  is  used  to  indicate  that  Ea{S)  is  the  contribution  to  the  interaction  energy 
E{S)  due  to  the  resolved  shear  stress  in  the  direction  of  the  Burgers  vector  r  =  [b  •  (a-^  •  n)]/6  on 
the  glide  plane  S. 

We  wish  to  determine  the  nodal  force  Gf  which  is  given  by  the  rate  of  change  of  Ea^S)  with 
respect  to  the  position  of  the  approximate  dislocation  configuration.  Unlike  the  stress  field  cr^, 
the  applied  stress  field  is  arbitrary  and  is  one  of  the  input  variables  of  the  model.  In  most 
cases,  it  is  necessary  to  evaluate  the  integral  in  (A.27)  numerically.  The  simplest  quadrature  rule 
is  obtained  by  assuming  that,  for  each  dislocation  segment  Cj,  r  is  a  function  of  position  along  the 
segment  and  depends  linearly  on  the  values  at  the  nodes.  In  other  words, 

r(£)  =  Tj-i  +  -(r,-  -  Tf-i)  0  <  £  <  li,  (A.28) 

where  t,_i  and  rj  are  the  values  of  r  at  the  nodes  a,_i  and  a;  of  the  segment  C,  of  length  1;,  and 
the  sense  of  the  arclength  i  is  from  node  a,_i  to  node  n,-. 

In  light  of  (A.28),  it  can  be  shown  that  the  rate  of  change  of  Ea{S)  witli  respect  to  to  a  change 
in  Of,  while  holding  all  other  nodes  fi.xcd,  is  given  by 

Gf  =  +  ^{di  +  di+i)Ti  +  ^d.+iri+ij,  (/1.29) 

where  d,  and  dj+i  are  the  horizontal  distances  between  the  nodes  fl,_i  and  «,  and  between  a,  and 
a, -4.1,  respectively. 
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FIGURE  CAPTIONS 


Fig.  1  A  threading  dislocation  gliding  down  the  a:i,X2-plane  between  an  embedded  strained  layer 
of  thickness  h  and  two  half-spaces.  The  threading  dislocation  extends  from  the  half-spaces 
through  the  layer,  and  it  then  continues  along  each  of  the  interfaces  as  a  mismatch  dislocation. 

Fig.  2  A  discrete  approximation  of  the  actual  configuration  of  Figure  1  consisting  of  TV  -f  2  linear 
dislocation  segments.  The  numerical  simulations  are  carried  out  in  the  xi,X2-plai^e,  the  glide 
plane  of  the  threading  dislocation. 

Fig.  3  The  simplest  of  discrete  approximations  to  the  configuration  of  Figure  1.  The  glide  plane  of 
the  dislocation,  the  ii,X2-pl3.ne,  is  inclined  at  an  angle  a  from  the  normal  to  the  bi-material 
interface. 

Fig.  4  Two  dislocation  segments  C,  and  C,,i  which  are  part  of  the  piecewise  configuration  of  Figure 
2  in  a  local,  oblique  x,2/-coordinate  system. 

Fig.  5  The  motion  in  the  xi,X2-pl^i^®  of  tho  threading  dislocation  at  an  early  non-steady  stage  of  its 
growth  history.  The  thickness  of  the  embedded  layer  is  1.1  times  the  critical  thickness. 

Fig.  6  The  long-term,  steady-state  motion  of  the  threading  dislocation  in  the  Xi,X2-plano.  The 
thickness  of  the  embedded  layer  is  1.1  times  the  critical  thickness. 

Fig.  7  Non-dimensional  values  of  the  nodal  driving  force  G,  for  various  nodes  in  the  layer  as  a  function 
of  time.  The  time  normalization  factor  to  is  equal  to  {b/ZVo)exY>{Qo/kT).  The  labels  on  the 
curves  correspond  to  xj  coordinates  of  the  nodes  in  terms  of  the  projected  hayer  thickness  hp. 
The  embedded  layer’s  thickness  is  1.1  times  the  critical  thickness. 

Fig.  8  Non-dimensional  values  of  the  nodal  driving  force  G,  for  various  nodes  in  the  layer  as  a  function 
of  nodal  position.  The  labels  on  the  curves  correspond  to  xi  coordinates  of  the  nodes  in  terms 
of  the  projected  layer  thickness  hp.  The  thickness  of  the  embedded  Layer  is  l.l  times  the  critical 
thickness. 

Fig.  9  The  motion  in  the  Xi,X2-plane  of  the  threading  dislocation  at  an  early  non-steady  stage  of  its 
growth  history.  The  thickness  of  the  embedded  layer  is  two  times  the  critical  thickness. 

Fig.  10  The  long-term,  steady-state  motion  of  the  threading  dislocation  in  the  Xi,X2-plane.  The 
thickness  of  the  embedded  layer  is  two  times  the  critical  thickness. 
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Fig.  11  Non-dimensional  values  of  the  nodal  driving  force  (j,  for  various  nodes  in  the  layer  as  a  function 
of  time.  The  time  normalization  factor  to  is  equal  to  {b/ZVo)exi>{Qo/kT).  The  labels  on  the 
curves  correspond  to  xi  coordinates  of  the  nodes  in  terms  of  the  projected  layer  thickness  hp. 
The  thickness  of  the  embedded  layer  is  two  times  the  critical  thickness. 

Fig.  12  Non-dimensional  values  of  the  nodal  driving  force  G,  for  various  nodes  in  the  layer  as  a  function 
of  nodal  position.  The  labels  on  the  curves  correspond  to  xi  coordinates  of  the  nodes  in  terms 
of  the  projected  layer  thickness  hp.  The  thickness  of  the  embedded  layer  is  two  times  the 
critical  thickness. 

Fig.  13  An  embedded  strained  layer  and  two  intersecting  glide  planes.  The  threading  dislocation  glides 
on  the  a:2 » ^a-pla^ne  in  the  positive  a;2-direction.  A  misfit  dislocation  lies  along  the  xi-axis. 

Fig.  14  Pyramids  of  the  close-packed  {111}  planes  in  a  cubic  crystal  showing  the  various  glide  planes. 
The  crystallographic  direction  normal  to  the  interface  is  [001].  The  Burgers  vector  of  the 
threading  dislocation  on  the  shaded  plane  coincides  with  either  of  the  edges  aq  or  bq,  while  the 
Burgers  vector  of  the  straight  dislocation  lying  along  ad  coincides  with  either  ap,  dp,  aq  or  dq. 

Fig.  15  The  motion  in  the  a:2,a:3-plane  of  a  threading  dislocation  that  was  started  at  X2  =  -1506.  A 
misfit  dislocation  is  positioned  at  *2  =  0  and  X3  =  /j/2cosa.  The  thickness  of  the  embedded 
layer  is  1.1  times  the  critical  thickness. 

Fig.  16  Non-dimensional  values  of  the  nodal  driving  force  G,  for  the  center  node  (xg  =  0)  and  the 
interfacial  node  {x'^/b  =  hp/2)  as  a  function  of  nodal  position.  The  dotted  curves  correspond 
to  the  driving  force  in  the  absence  of  the  blocking  misfit.  The  threading  dislocation  was  started 
at  X2  =  —1506. 

Fig.  17  A  close-up  of  the  motion  in  the  i2)®3-plane  of  a  threading  dislocation  that  was  started  at 
X2  =  —3006  in  the  vicinity  of  a  blocking  misfit.  The  misfit  dislocation  is  positioned  at  X2  =  0 
and  X3  =  /i/2cosa.  The  dotted  curves  correspond  to  the  motion  of  the  same  threading 
dislocation  in  the  absence  of  the  misfit.  The  thickness  of  the  embedded  layer  is  1.1  times  the 
critical  thickness. 

Fig.  18  Non-dimensional  values  of  the  nodal  driving  force  G',  for  the  center  node  (X3  =  0)  and  the 
interfacial  node  (X3/6  =  hp/2)  as  a  function  of  nodal  position.  The  dotted  curves  correspond 
to  the  driving  force  in  the  absence  of  the  blocking  misfit.  The  threading  dislocation  was  started 
at  *2  =  -3006. 
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Fig.  19  Non-dimensional  values  of  the  nodal  driving  force  G,  for  the  center  node  (.T3  =  0)  and  the 
interfacial  node  {x'^/b  =  hp/2)  as  a  function  of  nodal  position.  The  dotted  curves  correspond 
to  the  driving  force  in  the  absence  of  the  blocking  misfit.  The  threading  dislocation  was  started 
at  X2  =  -6006. 

Fig.  20  The  motion  in  the  a;2,a;3-plane  of  a  threading  dislocation  that  was  started  at  X2  =  -3006.  A 
misfit  dislocation  is  positioned  at  2:2  =  0  and  x'^  =  /i/2cosa.  The  thickness  of  the  embedded 
layer  is  two  times  the  critical  thickness. 

Fig.  21  Non-dimensional  values  of  the  nodal  driving  force  Gi  for  the  center  node  (2:3  =  0),  an  inter¬ 
mediate  node  (2:3/6  =  36p/8),  and  the  interfacial  node  (x'^/b  =  hp/2)  as  a  function  of  nodal 
position.  The  dotted  curves  correspond  to  the  driving  force  in  the  absence  of  the  blocking 
misfit.  The  threading  dislocation  was  started  at  X2  =  -3006. 

Fig.  22  Geometry  for  a  segment  Ci+i  interacting  with  a  second  disjoint  segment  Cm+i- 

Fig.  23  Geometry  for  a  segment  C,-  interacting  with  a  second  disjoint  segment  Cm+i. 

Fig.  24  Geometry  for  a  segment  Cj+i  interacting  with  a  second  adjacent  segment  Ci. 

Fig.  25  Geometry  for  a  segment  Cj+i  interacting  with  two  semi-infinite  segments  and 
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